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ABSTRACT 

We investigate why the spectral type of most cataclysmic variable (CV) secondaries 
is significantly later than that of a ZAMS star with the same mean density. Using 
improved stellar input physics, tested against observations of low-mass stars at the 
bottom of the main sequence, we calculate the secular evolution of CVs with low-mass 
donors. We consider sequences with different mass transfer rates and with a different 
degree of nuclear evolution of the donor prior to mass transfer. 

Systems near the upper edge of the gap (P 3 — 6 h) can be reproduced by models 
with a wide range of mass transfer rates from 1.5 x 10"^AfQyr~^ to 10~^MQyr~^. 
Evolutionary sequences with a small transfer rate and donors that are substantially 
evolved off the ZAMS (central hydrogen content 0.05 — 0.5) reproduce CVs with late 
spectral types above P ^ 6 h. Systems with the most discrepant (late) spectral type 
should have the smallest donor mass at any given P. 

Consistency with the period gap suggests that the mass transfer rate increases 
with decreasing donor mass for evolved sequences above the period gap. In this case, 
a single-parameter family of sequences with varying and increasing mass transfer 
rate reproduces the full range of observed spectral types. This would imply that CVs 
with such evolved secondaries dominate the CV population. 
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1 INTRODUCTION 

Over last two decades, numerous studies have considered 
the evolution of the low-mass secondaries in cataclysmic 
variable (CV) systems. The calculations aimed mainly at 
reproducing the observed distribution of orbital periods P 
of CVs, in particular the minimum period at 80 min and the 
dearth of systems in the 2-3 h period range (see e.g. King 
1988 for a review). Based on the most popular explanation 
for this period gap, the disrupted magnetic braking model 
(cf. Rappaport, Verbunt and Joss 1983; Spruit and Ritter 
1983), evolutionary sequences constructed with full stellar 
models or simplified bipolytrope models (cf. Hameury et al. 
1988; Hameury 1991; Kolb and Ritter 1992; and references 
therein) reproduce the broad features of the observed P dis- 
tribution. Despite the significant uncertainties in the under- 
lying stellar input physics (e.g. opacities, equation of state, 
etc.), the lack of a reliable description of angular momen- 
tum losses J that drive the mass transfer provided sufficient 
freedom that it was always possible to find a set of secondary 
star models which reasonably fit the period gap. 

Input physics like atmospheric opacities and the outer 



boundary condition are essential for a comparison with ob- 
servable quantities such as colours or spectral types. 

The spectral type and the orbital period [P oc 
{R^ /MY^'^; R and M are the secondary radius and mass) 
provide a unique set of constraints that test the structure 
both of the secondary's interior and its outermost layers. As 
we shall show below, these can be used to find constraints 
on the actual orbital braking strength J as well. 

Previously, an analysis of this kind was hampered by 
the lack of knowledge, both theoretical and observational, 
of stars in the low-mass range populated by CV secon- 
daries. Here we make use of the marked and independent 
progress in the study of low-mass stars over the past decade. 
A wealth of ground-based and spaced-based observations has 
significantly improved our knowledge of the photometric and 
spectroscopic properties of stars at the bottom of the main 
sequence. In parallel, recent theoretical work has demon- 
strated the necessity to use accurate internal physics and 
outer boundary conditions based on non-grey atmosphere 
models to correctly describe the structure and evolution of 
low-mass objects (Burrows et al. 1993; Barafle et al. 1995, 
1997, 1998; Chabrier and Baraffe 1997, 2000, and references 
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therein). These interior models, combined with recent, much 
improved models of cool atmospheres and synthetic spec- 
tra (see the review of AUard et al. 1997, and references 
therein) allow one to calculate self-consistent magnitudes 
and synthetic colours of low-mass stars. These can be com- 
pared directly to observed quantities, avoiding the use of 
uncertain empirical effective temperature T^s and bolomet- 
ric correction scales. The agreement between models and 
observations for (i) eclipsing binary systems (Chabrier and 
Baraffe 1995), (ii) mass-magnitude and mass-spectral type 
relations (Baraffe and Chabrier, 1996; Chabrier et al. 1996), 
(iii) colour-magnitude diagrams (Baraffe et al. 1997, 1998) 
and (iv) synthetic spectra (Leggett et al. 1996; AUard et 
al. 1997), demonstrate how reliable the theory of low-mass 
stars already is. 

Recently, Beuermann et al. (1998) used our prelimi- 
nary calculations based on these next -generation low-mass 
star models and compared calculated and observed spectral 
types (SpT) of CV secondaries. Beuermann et al. (1998) 
found that CVs with donors that are nuclearly evolved, 
with a substantial amount of central hydrogen depletion, 
can account for the late spectral types in long-period sys- 
tems (-P/h ^ 6), but not in those close to the upper edge of 
the period gap (3 ^ P/h 6). Higher mass transfer rates 
as usually assumed can account for the latter. 

The purpose of the present paper is to give a detailed 
account of the full set of our CV evolutionary calculations, 
a preliminary subset of which has been used by Beuermann 
et al. 1998. We discuss implications of the need for both 
high mass transfer rates and a significant fraction of systems 
with nuclearly evolved donors. We focus on CVs where the 
spectral type can be determined fairly accurately, i.e. on 
systems with orbital period above 3 h. In a previous paper 
(Kolb and Baraffe 1999), we considered aspects specific to 
the evolution of systems below the period gap (P < 2 h). 

In §2 we summarize the input physics and present evo- 
lutionary sequences with initial ZAMS donors. Nuclearly 
evolved sequences are presented in §3. A discussion and con- 
clusions follow in §4. 



2 EVOLUTIONARY MODELS WITH INITIAL 
ZAMS DONOR 

A brief description of the input physics for our stellar code 
can be found in Kolb and Baraffe (1999) and Baraffe et 
al. (1998). More details are given by Chabrier and Baraffe 
(1997) for the interior physics and by Hauschildt et al. (1999) 
for the atmosphere models. Here we just recall the main 
improvements compared to earlier models applied to CV 
donors: (i) the new models are based on the equation of state 
(EOS) of Saumon et al. (1995) which is specially devoted to 
low-mass stars and brown dwarfs, and to the description 
of strong non-ideal effects characterizing the interior stel- 
lar plasma; (ii) the outer boundary condition is based on 
non-grey atmosphere models (the use of a grey boundary 
condition overestimates the effective temperature Teff and 
luminosity L for a given mass M); (iii) mass-colour and 
mass-magnitude relations are derived self-consistently from 
the synthetic spectra of the same atmosphere models which 
provide the boundary condition. 

We use the empirical SpT - (I — K) relation estab- 



lished by Beuermann et al. (1998) to determine the spectral 
type SpT for a given model (M, Tes, L) from the calculated 
colour I — K. The accuracy and relevance of this conversion 
is discussed in Beuermann et al. (1998). 

Because of the uncertainty of current descriptions for 
orbital angular momentum losses by magnetic braking and 
the substantial computer time (cf. Kolb and Ritter, 1990; 
Hameury 1991) required for evolutionary sequences if mass 
transfer is treated explicitly and self-consistently, we re- 
strict our analysis to sequences calculated with constant 
mass transfer rate. We have checked that sequences with 
mass transfer driven by angular momentum losses J with 
either constant J, constant J / J, or J according to Verbunt 
& Zwaan (1981), have the same main properties as the se- 
quences with constant M on which we base our conclusions. 

In this section we consider sequences where mass trans- 
fer starts with donors of mass M2 = IMq that are ini- 
tially on the ZAMS, with solar composition (X = 0.70, 
Z=0.02). We define a "standard sequence" which reproduces 
the width and location of the period gap (2.1 - 3.2 h) in the 
framework of the disrupted magnetic braking model. This 
sequence is calculated with M = 1.5 x lO~^M0yr~^ un- 
til the donor becomes fully convective at M2 = 0.21 Mq 
and P = 3.2 h. At this point mass transfer ceases and the 
secondary shrinks back to its thermal equilibrium configu- 
ration in 2 — 3 X 10* years. Mass transfer resumes when the 
Roche lobe radius catches up with the donor radius at a 
period P = 2.1 h. The transfer rate in this second phase is 
5 X lO~^^M0yr~^, a value typical for systems driven by grav- 
itational wave emission (cf. Kolb and Baraffe 1999). This 
standard sequence is shown in Fig. ^ (thick solid line), to- 
gether with sequences calculated with higher transfer rates 
(3 X lO"^M0yr-\ lO"**M0yr~^ and lO"'^M0yr~^). 

Mass-losing stars deviate from thermal equilibrium. 
The deviation is large if the secondary's thermal timescale 
fKH ~ GkP / RL is long compared to the mass transfer 
timescale tu = M/{—M). As shown in Fig. mass trans- 
fer causes a contraction of the donor radius (e.g. Whyte 
and Eggleton 1980; Stehle et al. 1996) relative to the corre- 
sponding ZAMS radius for predominantly radiative objects 
M ^ O.6M0, where the radiative core exceeds 80% of the 
total mass (cf. Chabrier and Baraffe 1997, their Fig. 9). Pre- 
dominantly convective stars (M O.6M0) expand relative 
to the ZAMS. Although a large mass transfer rate causes 
a significant departure of R from its equilibrium value (see 
Fig. 1^), Fig. 1^ shows that T^s is rather insensitive to M 
(see also Singer et al. 1993; Kolb et al. 2000). 

The outer boundary condition of the stellar model is 
more important than M for the determination of T^s, as in- 
dicated by the dotted curve in Fig. |l] and Fig. |b. This curve 
corresponds to a sequence with M = 1.5 x lO~^M0yr~^, 
but calculated with the Eddington approximation, which 
implicitly assumes greyness and the diffusion approxima- 
tion for radiative transfer in the atmosphere. Chabrier and 
Baraffe (1997) have shown that such an approximation over- 
estimates Tcff for a given mass when molecules form in the 
atmosphere at T^s <, 4000 - 4500 K. Note that the Edding- 
ton approximation sequence reproduces the period gap as 
well as the corresponding non-grey sequence, with a width 
and location of 2.2 - 3 h. However, the Eddington sequence 
consistently gives earlier spectral types than observed, and 
dramatically fails to reproduce the location of objects with 
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Figure 1. Spectral type of CV secondaries versus orbital period, for evolutionary tracks with unevolved donor stars and different mass 
transfer rate. The observed data are taken from Beuermann et al. (1998), with slight updates. The short-dashed curve corresponds to 
solar metallicity stellar models on the ZAMS. The thick solid curve is the standard sequence with constant M = 1.5 X 10~^MQyr~^ 
above the period gap and M = 5 X 10~^^MQyr~^ below the gap. The thick long-dashed curve is with constant M = 3 X 10~®; the 
dash-dotted curve is with M = lO~*M0yr~^; the thin long-dashed curve is with M = lO~^M0yr~l. For all sequences, the initial 
secondary mass is M2 = 1 Mq. Sequences above the period gap terminate when the donor becomes fully convoctive. The dotted curve 
corresponds to a sequence with M = 1.5 X lO~®M0yr~^, but calculated with the Eddington approximation as outer boundary condition 
for the stellar models. 



P < 5 h in Fig. This example highlights the significant 
improvements of our new generation of low-mass star mod- 
els. 

Note that the period Ptum where the period derivative 
changes from negative to positive due to departure from 
thermal equilibrium increases with mass transfer rate. Re- 
markably, the standard sequence reaches period bounce at 
P = 3.2 h, the same period where the secondary star be- 
comes fully convective. Sequences with higher mass transfer 
rate bounce before the donor becomes fully convective. In 
other words, for mass transfer rates M > 1.5 x lO~^M0yr~^, 



the minimum period reached by CVs just corresponds to or 
exceeds the observed upper edge of the period gap at '--^ 3 h. 

As mentioned in Beuermann et al. (1998), a large 
spread of the secular mean mass transfer rates from 1.5 x 
lO"^M0yr"^ to ~ lO"*M0yr"^ could account for the ob- 
served range of late spectral type objects with P 5 h. It 
is well known that such a spread of M in systems above 
the period gap cannot be excluded by observations (cf. e.g. 
Warner 1995, his Fig. 9.8; Sproats et al. 1996). 

Although the sharp observed gap boundaries seem to 
imply uniform values of the secular mean M near the upper 
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Table 1. Donor mass M and mass Mcoro of its radiative core, at 
P = 3 h for evolved sequences with initial donor mass M2 = 1 
or 1.2 Mq, constant mass transfer M = 1.5 X lO^^MQyr"-*^ and 
different central H abundance Xc at the time mass transfer starts, 
to is the age of the donor when mass transfer turns on. 



M2/M0 




to (Gyr) 


A//M0 


A/corc/M0 


1 


0.55 


2.37 


0.191 


0.064 




0.16 


7.2 


0.209 


0.125 




0.05 


8.16 


0.199 


0.136 


1.2 


0.16 


3.75 


0.183 


0.120 




0.05 


4.18 


0.143 


0.127 



edge of the period gap (e.g. Ritter 1996), a spread of M can- 
not be dismissed unambiguously. The sharp upper edge is 
preserved if disrupted magnetic braking holds and M adopts 
only values > 1.5 x lO"^M0yr"\ With increasing M the 
critical mass where the secondary becomes fully convective 
becomes smaller (cf. Fig. ^). Thus, if the system detaches 
at this point, by analogy with the disrupted magnetic brak- 
ing model, mass transfer resumes after the secondary has 
reached thermal equilibrium at a period shorter than 2 h. 
As an example, for the sequence with M = lO~*M0yr'^ 
the donor becomes fully convective at M2 ~ O.I5M0, and 
the system reappears as CV at P ~ 1.7 h. Test calculations 
with a population synthesis code (e.g. Kolb 1993) show that 
such a spread of systems reappearing below the gap still 
gives a reasonably sharp lower edge of the period gap, with 
only a mild over-accumulation of systems in a period inter- 
val 1.7 — 2.1 h. This is perfectly consistent with the observed 
distribution which does show a slight, although statistically 
probably not significant, accumulation of systems there. 

Despite this consistency, accepting the M spread is not 
an attractive solution to the late SpT problem. It is hard to 
see why a presumably global angular momentum loss mech- 
anism should drive such different transfer rates in otherwise 
similar systems. 

Finally we note that there are no CV secondaries above 
the ZAMS line in the P— SpT diagram for periods longer 
than ~ 6 h. This suggests that in these systems the secular 
mean M must be smaller than lO~*M0yr~^. But whatever 
the value of M, sequences with unevolved donors never drop 
below the ZAMS line for P ^ 6 h. 



3 EVOLVED SEQUENCES 

We now consider sequences with initial donors which have 
evolved off the ZAMS prior to mass transfer. When mass 
transfer turns on, the star is old enough (cf. Table 1) to have 
burned a substantial amount of hydrogen in the centre. Only 
stars with mass AI > O.35M0 , which develop a radiative core 
(cf. Chabrier and Baraffe 1997), start to deplete central H 
within a Hubble time. But only if M > O.8M0 is more than 
50% of the initial H depleted in less than 10 Gyr. Therefore, 
this donor type is restricted to systems that form with a 
fairly massive secondary ^ IMq. 

Figure ^ depicts sequences with such nuclearly evolved 
donors — hereafter "evolved sequences" — in the P-SpT 
diagram, for different central H mass fractions Xc at the time 
mass transfer starts. The initial secondary mass is AI2 = 



IM0 or M2 = I.2M0 and mass transfer rate is fixed to 
1.5 X lO~®M0yr~^. Table 2 gives for each of these sequences 
the secondary mass A4, spectral type, T^h and radius P as a 
function of P, for 3 < P/h < 10. The evolutionary track of 
ZAMS donor sequences in the P-SpT diagram is insensitive 
to the initial donor mass, while for evolved sequences with 
the same initial Xc Teff decreases (the spectral type becomes 
later) at a given P with increasing initial donor mass. This 
is illustrated in Fig. ^ where we plot sequences starting with 
M2 = IM0 (long-dashed curve) and M2 = I.2M0 (dash- 
dotted curve), both for Xc = 0.05. 

As already emphasized by Beuermann et al. (1998), 
evolved sequences result in later spectral types for a given 
P ^ 5 h than the standard sequence and can explain amaz- 
ingly well the spread of data observed above 6 h. The reason 
for this is illustrated in Figures ^ and |^, which display sev- 
eral diagnostic quantities along the sequences as a function 
of secondary mass (Fig. |^ and orbital period (Fig. |^ . 

When mass transfer turns on, more evolved donors with 
smaller Xc have larger P, Tcft and Tc than less evolved 
donors. This follows from the well known properties of the 
If burning phase of solar-type stars with a central radiative 
core: as hydrogen is depleted in the centre, the increase of the 
central molecular weight fi yields an increase of the central 
temperature and thus L, causing an expansion of the star. 
Hence the main-sequence evolution of a star with constant 
mass proceeds towards larger L and Tcb . The same effect ex- 
ists for mass-losing main-sequence stars. Consequently, for 
a given mass, R is larger for more evolved sequences, and P 
is longer, while for a given P, the secondary mass and Tcfi 
is smaller, hence the spectral type is later (cf. Fig. As a 
consequence, the donor mass at any given P is smallest for 
the most evolved sequence, cf. Table 2. 

Although the evolved sequences shown in Fig ^ provide 
an attractive explanation for the observed data above 5-6 
h, they are obviously in confiict with the standard expla- 
nation of the period gap. Indeed, with the adopted value 
M — 1.5 X lO~^M0yr~^ we find that these sequences en- 
ter the period gap before the donor becomes fully convec- 
tive or before the system bounces due to departure from 
thermal equilibrium. Table 1 gives the mass of the sec- 
ondary and the mass Mcore of its radiative core at P = 
3 h for the different sequences. Period bounce occurs for the 
(M2 = IM0, Xc = 0.16) sequence at P = 2.60 h, corre- 
sponding to M = O.IM0 and Mcoro = O.O4M0. The donor 
becomes fully convective only at M ~ O.O5M0 (P = 2.9 h). 
For more evolved sequences, period bounce would occur at 
even smaller periods. 

This difference in the size of the radiative core is caused 
by the small central H abundance. The smaller Xc and the 
higher central temperature of evolved donors compared to 
unevolved donors (cf. Figs. ^ and ^) imply lower radiative 
opacities in the central region, and consequently favour ra- 
diative transport. As a result, the inward progression of the 
bottom of the donor's convective envelope as mass decreases 
proceeds at a slower rate. Hence in evolved sequences the 
donor has a larger radiative core Mcoro than a donor with 
the same mass in an unevolved sequence (see Fig. ^). In 
terms of P, Fig. ^ shows that Mcore is larger for the stan- 
dard sequence at any given P ^ 3.5 h. The situation then 
reverses at the upper edge of the gap, which corresponds to 
Mcore ~ for the standard sequence. The effect of a lower 
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Figure 2. (a) Mass - radius relation for ZAMS stars and for unevolved CV donors in sequences with constant mass transfer rate. 
Linestyle as in Fig. 1. (b) Mass — T^g relation for the same models as in (a). 



central H abundance described above is similar to that re- 
ported by Laughlin et al. (1997) who followed the central 
H burning phase of very low-mass stars with M < O.25M0 
(with constant mass). They note that stars which are fully 
convective on the ZAMS develop a radiative core towards 
the end of the main sequence because of the depletion of H 
and subsequent lowering of the central radiative opacities. 

The P-SpT diagram suggests that at long orbital pe- 
riods a fairly large fraction of systems has a significantly 
evolved donor. If these carried on evolving with a trans- 
fer rate ~ 1.5 x lO~^M0yr~^ they would clearly overpop- 
ulate the period gap regime. This is made worse by the 
fact that some of these sequences reach their minimum pe- 
riod in the period gap (Fig. at period bounce P = 0, 
and the detection probabihty d oc 1/P increases sharply. 
We find a similar behaviour if mass transfer is treated self- 



consistently, with an angular momentum loss rate J — const, 
or J (X J, calibrated such that the corresponding un- 
evolved sequence reproduces the observed period gap (giving 
M = 1.5 X IQ-'^Meyr-i at P ~ 3 h). 

An obvious way to prevent evolved sequences from 
evolving into the period gap region is to increase the mass 
transfer rate. This increases the deviation from thermal equi- 
librium and leads to period bounce above the upper edge 
of the observed gap. Test calculations with constant mass 
transfer rate show that period bounce occurs at P > 3 h 
with M > 3 X 10~^ and 5 x 10~^ MQyr~^ for sequences 
with respectively Xc = 0.16 and 0.05, for initial donor mass 
IMq. The same behavior is found for donors with higher 
initial mass. In other words, for sequences with smaller Xc a 
higher M near the upper edge of the period gap is required 
than for less evolved donors, consistent with the fact that 
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Figure 3. As Fig. 1, but for sequences with different initial central H abundance Xc at the time mass transfer starts. The mass transfer 
rate is M = 1.5 X 10~^MQyr~^, the initial secondary mass M2 = 1 and 1.2 M©. Thick solid line: standard sequence with M2 = IMq 
and Xc = 0.7; dotted curve: M2 = IMq and Xc = 0.55; dashed curve: M2 = IMq and Xc =il.l6; long— dashed curve; M2 = IMq and 
Xc = 0.05; dash-dotted curve: M2 = I.2M0 and Xc = 0.05. The data are the same as in Fig. hi 



the thermal time is shorter if Xc is small. These evolved se- 
quences with high transfer rates {M > 5 x lO~^M0yr~^) 
also match the observed spectral types in the region P = 
3-6 h, just as high M sequences with ZAMS donors do (see 
§2 and Fig. §). 

Note however, that for long periods P ^ 5 h a higher 
M has the same effect as for the unevolved sequences dis- 
cussed in §2 (Fig. The effective temperature is larger and 
therefore the spectral type earlier than for a sequence with 
smaller M. Thus in order to account for the full range of 
rather cool spectral types at long periods, M cannot always 
be high along all evolved sequences. 

To summarize: evolved sequences require high M (^ 
5 X 10~®) near the upper gap region to avoid crossing the 



gap, but a low mass transfer rate [M ^ 5 x 10"'') for P ^ 6 
h to explain the late spectral type CVs in this region. This 
suggests that the mean mass transfer rate increases during 
the secular evolution of nuclear evolved donors. 



4 DISCUSSION AND CONCLUSIONS 

We performed calculations of the long-term evolution of 
CVs with up-to-date stellar models that fit observed prop- 
erties of single low-mass stars exceptionally well. Our focus 
was on the orbital period- spectral type (P-SpT) diagram, 
and the observed location of CV secondaries which populate 



© 0000 RAS, MNRAS 000, 000-000 



On the late spectral types of cataclysmic variable secondaries 7 




0,2 0.4 0.6 0.8 1 1.2 

M/M, 

Figure 4. Variation of several stellar quantities with mass 
for sequences with constant mass transfer rate M = 1.5 x 
lO~^M0yr~i, initial masses M2 = 1 and 1.2 Mq and differ- 
ent central H abundances Xc at the time mass transfer starts. 
Solid curve: standard sequence with M2 = IMq and Xc = 0.7; 
dashed curve: M2 = IMq and Xc = 0.16; long-dashed curve: 
M2 = IMq and Xc = 0.05; dash-dotted curve: M2 = 1.2Mo 
and Xc = 0.05. Mcorc is the mass of the radiative core. Tc is the 
central temperature (in K). 

a band with spectral types significantly later than ZAMS 
stars at the same orbital period. 

Our calculations tested if the observed spectral types 
can be reproduced by varying two main parameters, the 
mass transfer rate M (assumed constant), and the degree 
of nuclear evolution of the secondary before mass transfer 
starts (measured in terms of the initial central H fraction 
Xc). Only CVs forming with massive (^ 1 — I.2M0) donors, 
i.e. at long periods, can have donors where H is significantly 
depleted in the centre. 




5 10 15 

P(h) 

Figure 5. Same as Fig. 4, but as a function of orbital period. 



We summarize our results as follows: (i) In the frame- 
work of the standard discontinuous orbital braking model 

for the period gap, the observed gap between 2-3 h is well 
reproduced by sequences starting from ZAMS donors and 
proceeding at a mass transfer rate M ~ 1 — 2 x lO^'^Moyr"^ 
near the upper edge of the gap. (ii) Higher transfer rates 
(M ^ 5 X 10~^Moyr~^) than assumed in (i) near the up- 
per edge of the gap arc required to explain the full range 
of late spectral types of secondaries in CVs with periods 
P = 3 — 6 h. This is true whatever the evolutionary state of 
the donor at turn-on of mass transfer, (iii) A family of evo- 
lutionary sequences starting mass transfer from an evolved 
donor with varying Xc < 0.5 and low mass transfer rate 
(M ^ 5 X IQ-^Moyr-i) covers the observed locations of 
CVs in the P— SpT diagram for P ^ 6 h. (iv) For these se- 
quences, higher mass transfer rates than in (i) are required 
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Table 2. Characteristic quantities as a function of P (for 3 < 
P/h < 10) for sequences with initial donor mass M2 = 1 or 
1.2 M0, constant mass transfer rate M = 1.5 X lO~^A/0yr~-^ 
and different central H abundance Xc at the time mass transfer 
starts. 



Mo 


X 


p 


M 


SpT 


T a 
^ eft 


R 


iV10 




h 






K 


J 1,0 


1.00 


0.700 


9.0 


1.000 


G5-G6 


5360. 


1.0441 






8.0 


0.992 


G6 


5289. 


0.9454 






7.0 


0.870 


KO-Kl 


4870. 


0.8205 






6.0 


0.751 


K4-K5 


4467. 


0.6972 






5.0 


0.616 


MO-MI 


3979. 


0.5758 






4.0 


0.475 


M2-M3 


3648. 


0.4574 






3.5 


0.377 


M3 


3521. 


KJ .OOOU 


1.00 


0.550 


9.0 


0.976 


G6 


5276. 


9933 






8.0 


0.912 


G8 


5105. 


0.9111 






7.0 


0.795 


K1-K2 


4713. 


0.7784 






6.0 


0.684 


K6 


4293. 


0.6666 






5.0 


0.571 


Ml 


3905. 


0.5626 






4.0 


0.437 


M2-M3 


3632. 


0.4433 






3.0 


0.191 


M4 


3324. 


0.2764 


1.00 


0.160 


10.0 


0.891 


G7 


5201. 


1.0617 






9.0 


0.815 


G9-K0 


4961. 


0.9571 






8.0 
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near 3 h, otherwise they would evolve into the period gap 
and predict too early spectral types at shorter P. 

If M is sufficiently large, the sequences reach their min- 
imum period above the upper edge of the period gap. The 
donor becomes fully convective at a mass smaller than the 
canonical O.2IM0 of the standard (unevolved) sequence, but 
at a period longer than 3 h. If the system detaches at this 
point and the donor can re-establish thermal equilibrium, 
then mass transfer would resume at a period well below 2 h. 
This guarantees consistency with the standard discontinu- 
ous orbital braking model. 

Points (ii) to (iv) show that if the mass transfer rate 
increases significantly along the secular evolution of evolved 



donors, these sequences can explain the observed scatter in 
the SpT diagram both above 6 h, through the effect a 
lower Xc has on the evolutionary track (see §3 and Fig. 
and below 6 h, through the effect of a higher AI (see §2 and 
Fig. 0). 

A large range of the secular mean mass transfer rate, as 
suggested by (ii) , could be reconciled with the sharp bound- 
aries of the observed period gap. In the standard model 
these are a result of the dominance of sequences with the 
canonical value 1 — 2 x lO~^M0yr~^ at periods close to 3 h. 
The apparent scatter of data in the P-SpT diagram does 
not support this dominance. Reasonably sharp boundaries 
could be maintained if the canonical value is the minimum 
M value of the range of transfer rates. The observed gap 
edges would then be defined by this minimum mass trans- 
fer rate sequence. Sequences with higher M would bounce 
before they reach the upper edge of the gap. They detach 
at longer periods and reattach at shorter periods than the 
standard sequence. 

A more serious problem is that the large M range for 
otherwise similar system parameters, as suggested by (ii) for 
unevolved sequences, leaves the nature of the main control 
parameter determining the strength of orbital angular mo- 
mentum losses completely undetermined. Our experiments 
show that Xc could be this parameter. 

The full observed range of data in the P-SpT diagram 
could be explained if the orbital braking strength increases 
both with decreasing P and decreasing Xc . This ensures that 
period bounce prevents evolved sequences from overpopulat- 
ing the period gap, while they still pass through rather late 
spectral types at long periods. The resulting J law must give 
the canonical value of AI at 3 h for Xc — 0.7. A large M 
range for unevolved systems is not required. 

With standard magnetic braking laws (e.g. by Verbunt 
& Zwaan 1981, Mestel & Spruit 1987) the transfer rate 
usually increases with period, and is smaller for sequences 
with more evolved donors (see e.g. Pylyser & Savonije 1989, 
Singer et al. 1993, Ritter 1994), i.e. has just the opposite 
differential behaviour than the one we propose. Given the 
significant uncertainty in these magnetic braking models 
and our poor knowledge of the underlying stellar magnetic 
dynamo, the suggestion that observations favour a non- 
standard law does not appear unrealistic. Further considera- 
tions of the theory of magnetic braking to assess our finding 
are clearly beyond the scope of this paper. We note, however, 
that there are magnetic braking scenarios in the literature 
that lead to an increase of M with progressing evolution. 
Kolb & Ritter (1992) obtained this for a Verbunt & Zwaan 
(1981)-type law when only the donor's convective envelope 
is coupled to the orbital motion. Zangrilli et al. (1997) found 
this property for their boundary-layer dynamo. 

At this point we re-emphasize that the input physics on 
which our stellar models are based provides an excellent de- 
scription of observed properties of isolated low-mass stars. 
Unless there are effects peculiar to CV secondaries that are 
not taken into account by the models (e.g., effects due to 
the rapid rotation, or the irradiation by the white dwarf 
and accretion disc) , one is forced to accept that a large frac- 
tion of CVs have a nuclearly evolved donor star. This is 
in stark contrast to predictions of standard models for the 
formation of CVs, where most CVs form with a secondary 
which is too low-mass (^ O.6M0) to be evolved (Politano 
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1996, King ct al. 1994). But this predominance of uncvolvod 
donors is mainly due to the neglect of systems forming with 
a secondary that is significantly more massive than the white 
dwarf (Ritter 2000). These undergo thermal-timescale mass 
transfer at a rate ^ lO^^Meyr^^ and arc usually associ- 
ated with supersoft sources (e.g. di Stefano and Rappaport 
1994). It is perfectly possible that these systems reappear 
as standard CVs once the mass ratio {q = donor mass/WD 
mass) is sufficiently small. In descendants of this thermal- 
timescale evolution any degree of nuclear evolution prior to 
mass transfer is possible. De Kool (1992) finds that the im- 
pact of these survivor systems depends on the assumed ini- 
tial mass ratio distribution which ultimately determines the 
secondary mass distribution in post common envelope bina- 
ries. For de Kool's model 1 {dN ocdg) survivor systems com- 
pletely dominate the CV population. A full appraisal of the 
viability of the apparent predominance of nuclear evolved 
systems is beyond the scope of our study. This requires a 
full population synthesis with self-consistent treatment of 
Xc-dependent angular momentum losses. 

Finally, the accurate determination of donor masses 
could provide an observational test of our hypothesis that a 
significant fraction of CV donors is nuclear ly evolved. Ob- 
servational errors arc still too largo (e.g. Smith & Dhillon 
1998) for this test. Systems with the most discrepant (late) 
spectral type should have the smallest mass at any given P. 
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